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Approximating  tail  probabilities  of  noncentral  distributions 


Suojin  Wang  and  H.  L.  Gray 
Department  of  Statistical  Science 
Southern  Methodist  University 


Abstract 

Gray  and  Wang’s  (1991)  general  method  for  approximating  tail  probabilities  is  applied  to  the 
cases  of  noncentral  xJi  F  and  t  distributions.  The  validity  of  such  applications  is  established.  The 
resulting  approximations  are  easy  to  compute.  Numerical  results  show  the  great  accuracy  of  the 
approximations  for  all  three  most  commonly  used  noncentral  distributions. 

Keywords:  Gp  ^-transform,  noncentral  \2,  noncentral  F,  noncentral  t,  tail  probability. 


1.  Introduction 

Calculating  the  tail  probabilities  of  noncentral  distributions  such  as  noncentral  \2,  F  and  t  is  a 
major  step  in  many  statistical  applications.  However,  these  distributions  do  not  have  a  simple  closed 
form  for  the  cumulative  distribution  functions  (CDFs)  which  may  be  expressed  in  the  form  of  infinite 
series  only.  Many  methods  have  been  proposed  for  approximating  the  CDFs  of  these  noncentral 
distributions  by  truncating  the  infinite  series  or  using  central  distributions  as  approximations.  These 
methods  are  designed  for  each  specific  distribution. 

Taking  a  totally  different  approach  from  the  viewpoint  of  the  Generalized  Jackknife,  Gray  and 
Wang  (1991)  have  introduced  a  general  method  referred  to  as  the  Gq  ^-transformation  for  finding 
functions  which  are  easy  to  evaluate  and  give  very  good  approximations  to  tail  probabilities.  The 
transformation  is  based  only  on  a  general  class  of  differential  equations  that  include  the  specified 
density  function  in  the  solution  set.  As  a  result  this  method  applies  to  a  broad  class  of  distributions. 

Since  the  densities  of  noncentral  distributions  are  themselves  in  the  form  of  infinite  series,  it  is  far 
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from  obvious  how  the  conditions  for  the  Gq  -transformation  are  met  in  these  cases.  In  this  note  we 
show  that  the  three  most  commonly  used  noncentral  distributions  (noncentral  xJ,  F  and  t)  all  satisfy 

(m ) 

the  conditions  so  that  the  Gq  -transformation  provides  a  good  alternative  in  approximating  the  tail 
probabilties  of  these  distributions.  Its  high  degree  of  accuracy  even  in  the  extreme  tails  is  maintained, 

as  will  be  shown  in  numerical  examples. 

(m) 

We  first  briefly  review  the  G^  -transformation.  Let  f  be  a  density  function  and  assume  that  we 
wish  to  approximate  the  tail  probability 

,°o 

G(x)  =  j  f(t)dt.  (1) 


f|,  00  ab  j 

Lk(x)  =  x k  z  ~r  ’ 

i=0  x1 

where  q  ^  0  and  is  an  integer  with  fj.  <  k.  Suppose  that  f(x)  satisfies  the  following  differential 


equation 


b m ( x )  /"'’(x)  +  Um_,(x)  l\x)  +  .  .  .  +  U j (x)  f^x)  -  f(x)  =  0  , 


for  some  collection  of  l'k's.  where  m  is  assumed  to  be  the  smallest  possible  integer  such  that  (2)  holds. 

( ni ) 

Using  the  idea  of  the  Generalized  Jackknife,  Gray  and  Wang  (1991)  define  the  Gn  -transformation 
approximating  G(x)  in  (1)  as 


G„mV(x);  ajj(x)]  = 


0 

f(x) 

f  (x) 

an(x) 

a12(x) 

ai,N+i(x) 

aN,j(x) 

aN,2<X) 

aN,N+t(x) 

1 

0 

0 

au(x) 

an(x) 

ai.N+,<x) 

aN,i -x^ 

aN,2(x) 

aN,N  +  ,W 

I H 


where 


As  shown  in  Gray  and  Wang  (1991),  most  commonly  used  simple  continuous  distributions 
including  those  in  the  Pearson  family  satisfy  (2)  for  m  =  1.  Furthermore,  it  has  been  shown  that 
under  mild  regularity  conditions, 


GnmV(x);  ajj(x)]  =  G(x)  . 

Gray  and  Wang  (1991)  have  examplified  the  extraordinary  accuracy  of  the  G 
small  as  2  or  3. 


(4) 


(m) 

n 


-transform  for  n  as 


2.  Noncentral  \2  distribution 

( IT1 ) 

The  first  application  of  the  Gn  -transformation  we  consider  is  to  the  noncentral  chi-squared 
distribution  \ ^ 2( A),  where  k  is  the  degree  of  freedom  and  A  is  the  noncentrality  parameter.  Cox  and 
Reid  (1987)  obtain  an  approximation  to  the  distribution  of  x^(^)  by  adding  a  small  perturbation  to 
the  corresponding  central  random  variable  when  A  is  small  relative  to  k.  Cohen  (1988)  and  Poeten 
(1989)  have  developed  algorithms  for  the  distribution  of  x^(^)-  These  algorithms  require  evaluating 
either  the  central  x£  distribution  function  or  the  distribution  functions  of  x^A)  for  the  lowest  degrees 
of  freedom  i  =  1,  2  and  3. 

The  density  of  X^(^)  can  be  expressed  as  (Johnson  and  Kotz,  1970,  p.  132) 


-0+10/2  OC  (1) 

«»>  =  e-  —  z  Ul 


X\T  xk/2  +  r-i 


2k/2  r=0 


r!  T(k/2  +  r) 


,  x  >  0. 


(5) 
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Now  note  that 


f(x)  =  e  gk(x)  0F!  !j;  , 


where  gk(x)  =  e  ’/(iXjj)  2^*)  is  the  density  of  , 


°Fl(  1  b!  Z)  “r5o  r!(b)r 


is  a  hypergeometric  function  (Rainville,  1960,  p.  74)  and 


(b)r  = 


1  ,  r  =  0  . 

b(b+l)  .  .  .  (b+r-1).  r  =  1,  2, - 


Let  u  =  0 F i ( — :  b;  z).  It  has  been  shown  (Rainville,  1960.  p.  109)  that 


zpi  +  b^M 

dz2  dz 


Let  b  =  ~  and  z  =  It  follows  from  (7)  that 


4x  d2u  i  2k  du  _  n 
A  dx2  +  A  dx  ~  U  ‘ 


Substituting  u  =  e  '  f(x)/g^(x)  into  (8)  and  after  some  algebra  we  have 

X  f"(x)  +  (x  +  (2-|))  f'(x)  +  *-±4~k~A  f(x)  =  0  ,  (9) 

where  we  used  the  fact  that  g|,(x)  =  ^  — ^  +  ()j  —  l)/xj  g^(x).  We  have  therefore  shown  that  f(x) 

(^) 

satisfies  (2)  with  m  =  2  and  £,  =  £2  =  0.  Hence  the  application  of  the  -transform,  denoted  by 
Gn  ^[f(x);  a—(x )],  is  straightforward  as  long  as  /  \x),  r  =  0,  1,  .  .  .  ,  2n+l,  are  obtained. 

Taking  r-th  derivative  of  (9)  and  after  some  algebra  we  have  the  recursive  formula 
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/r+l)(x)  =  |g  -J— x)  +  (fet*=i  -1-,) !»'»(,)  -  5  (<'",|(X)|/X  , 


(10) 


for  r  =  1,  2,  .  .  .  .  From  (9)  it  is  easily  seen  that  equation  (10)  is  also  valid  for  r  =  0  if  we  define 
/  \x)  =  0.  Thus  we  need  only  to  evaluate  f(x)  and  f^x).  But  from  (6)  we  have 

*V)  =  (-  5  +  *£2)  f(x)  +  ^  e“A/2gk(x)  0Fj(  — ;  if)  .  (11) 

By  (6)  and  (11),  it  is  sufficient  to  compute  a  hypergeometric  function  of  the  form 


oc  r  oo 

qFi(-;  b;  z)  =  £  riTry  =  £  cr, 
r=0  r-*b)r  r=0 


where  c0  =  1.  cr  =  — — — — —  cr_,,  r  >  1.  The  series  in  (12)  converges  as  quicklv  as  Y~  -^=—5 
r(b+r~l)  1  1  _  0  yr'j  (r!)2 

asymptotically.  Note  that  Levine  and  Sidi’s  (1981)  d-transformation  is  useful  to  accelerate  the 

convergence  of  the  series  in  (12). 

(2) 

To  compute  Gn  [f(x);  a-j(x)]  in  (3)  a  subroutine  for  computing  the  determinant  of  a  matrix  is 
normally  required.  But  when  n  =  1  it  is  easily  derived  that  in  the  noncentral  \2  case. 


Gl2)[f(x);  aij(x)]  =  ,  (13) 

where 

A(x)  =  “*+«)“  -3)]  f*(x) 

+  [|(-5x+A+5k~20)-(|  -2)(|  -3)Jf2(x)f/(x) 

+  [_2x2+(k-4)x]  f(x)^fr(x)^-x2(f'(x))3  , 
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B(x)  =  [s(-x+k+A-4)-^+2]  (f'(x))3+[f(-x+  ^+A-6)  +  |(k+A-4)(2-k)]f(x)f'(x) 

-  j^(  — x+k+A— 4)2  f2(x). 

(2)  (2) 

Table  1  provides  some  selected  values  of  Gq  .  Rapid  convergence  of  Gq  as  n  increases  is 
evident;  great  accuracy  is  obtained  for  n  as  small  as  2.  It  is  seen  that  the  approximation  is  particularly 
useful  for  extreme  tails.  We  used  a  finite  sum  to  approximate  the  series  in  (12),  and  therefore  f(x)  and 
f'(x).  In  Table  1,  15-40  terms  were  needed  depending  on  different  parameters  to  get  at  least  12 
significant  digits  correct  for  f(x)  and  f^x).  Applying  Levine  and  Sidi's  (1981)  acceleration  method 
improves  this  by  at  least  3  more  significant  digits.  However  the  quick  convergence  in  (12)  makes  it 
generally  unnecessary  to  use  Levine  and  Sidi's  method.  In  this  and  the  following  applications,  we  used 
the  IMSL  subroutine  LINV3F  to  calculate  determinants.  Any  other  efficient  algorithms  can  be  used  for 
this  purpose. 


3.  Noncentral  F  distribution 
( rn ) 

We  now  consider  the  application  of  the  Gn  -transform  to  noncentral  F  distributions.  Despite 
the  different  structures  of  the  noncentral  F  and  noncentral  \2  distributions,  a  similar  procedure  can  be 
developed.  Like  in  the  noncentral  \ 2  case  our  approach  is  different  from  those  usual  ones  which  use 
central  F  distributions  as  approximations;  see  Hirotsu  (1979).  The  density  of  a  noncentra)  F 
distribution  with  degrees  of  freedom  k  and  p  and  noncentrality  parameter  A  (denoted  by  F(  (A))  is 
(Johnson  and  Kotz,  1970,  p.  191) 


--V2/k>/2  k/2  -  1 


g(y)  = 


e  (§)“'*  y“'‘  ‘»ri^  +  r 


)rn 


fill  +  k,)(k+P>/!  ,!r,5  +  '>  V>+P>'. 


=  e 


r<2 

-A/2 


h  (y)  F  r ttE.  k.  a  _ky_ 

hk,p(y)  2  ’  2’  2  p+ky 


(14) 
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where 


r(^)  (&)k/V/,-‘ 

hk,p(y)  =  r(k) r(p)  {i+lrp+ p)/. 


is  the  density  for  &  central  F  distribution  with  degrees  of  freedom  k  and  p,  and 


‘F^a;  b:  z)  ~  rJQ  r!(b)r  Z' 


is  a  hypergeometric  function.  Since 


where  x  =  ^  y, 


/OO  ,OC  /-CO 

g(t)dt  =  J  g(^  s)  2  ds  =  J  f(s)  ds  , 
y  x  x 

*>- *)  I -•"VX(M  ■'■(*?= 


we  need  only  to  consider  the  transformed  density  f(x).  Let  u  =  jF,(a;  b;  z).  Using  the  fact  that 


(Rainville,  1960.  p.  124) 


z£?+(b-z)t  -au  =  ° 


and  after  lengthy  algebra  similar  to  those  in  Section  2  we  can  obtain  that 


«"<«>  -  ^  (1  -  ri)+ 


i-rb)+^  'w-  <“> 


Comparing  (18)  with  (2),  it  is  readily  seen  that  m  =  2,  lx  =  1  and  f2  =  2. 

A  recursive  formula  for  higher  order  derivatives  of  f(x)  can  be  easily  obtained  by  taking 
derivatives  on  both  sides  of  (18)  and  by  the  fact  that 


(h,(x)  h2(x))(0  =  £  (')  hij)(x)  h(2‘  J)(X),  i  =  1,  2 - 

j=0  J 
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It  is  therefore  only  necessary  to  evaluate  f(x)  and  f^x).  But 


f>,x)  _  fk-2  k+P~]  !■(,),  *  A/’(t+p)A 

W  L*  2(l+x)J  f(  1  +  2tr(p/2)(l+%,,k+p)/ 


„  /"k+p+2.  k+2. 

2  +  2  5  ’  2  ’ 


Ax  ) 
2(l+x)  r 


so  that  the  essential  part  is  to  compute  ^(a;  b;  2).  Similar  to  (12), 


iFj(a;  b;  z)  =  £  dr  , 
r=0 


(a+r  1  )z  00  r 

where  d0  =  1,  dr  =  — —  dr_!  ,  r  =  1,  2,  .  .  .  .  Therefore  ,F](a;  b;  z)  converges  as  fast  as  Z! 

asymptotically. 

(2) 

Selected  values  of  the  Gn  -transform  (n  =  1,  2,  3)  are  given  in  Table  2.  It  is  clearly  seen  that  n 

=  2  or  3  is  generally  sufficient  in  most  applications.  In  Table  2,  up  to  50  terms  in  (20)  were  used  to 

approximate  f(x)  in  (16)  and  f'(x)  in  (19)  to  at  least  12  significant  digits.  The  last  column  is  the 

Edgeworth-series  approximation  by  Mudhoikar,  Chaubey  and  Lin  (1976).  Table  1  shows  that  for  n  as 

(2) 

small  as  3  the  general  method  Gn  is  more  accurate  in  the  tail  than  the  Edgeworth  series 
approximation  designed  only  for  the  noncentral  F. 


4.  Noncentral  t  distribution 

This  section  concerns  with  noncentrai  t  distributions.  Kraemer  and  Paik  (1979)  proposed 
approximations  based  on  central  t  distributions  assuming  that  the  noncentrality  parameter  is  small 
relative  to  the  degree  of  freedom.  As  we  have  mentioned  earlier,  our  method  does  not  require  such 
assumption  on  the  noncentrality  parameter.  The  density  of  a  noncentral  t  distribution  with  k  degrees 
of  freedom  and  noncentrality  parameter  A  (denoted  by  tj^A))  is  (Johnson  and  Kotz,  1970,  p.  205) 


g(y)  = 


A2/2/  k  x(k  +  >)/2 

P^k  +  r  -f  1  ^ 

(&  A  y  V 

U+y2' 

£ r( ;  } 
r=0 

>f^k  f(jj) 

Vjk+y V 
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Let  x  =  1  +  y J/k.  It  follows  that 


/OO  .00  /  r  oo 

g(t)dt  =  j  g({k(s-l)}/) - t-_^d s  =  y  f(s)ds,  (22) 

y  x  2{k(8—  1)}  7  x 


where 


f(x)  =  g({k(x-l)}'/2) - - - - 


2{k(x  — I)} 


I)}1/2 


—  A2/2 


r(|')(x-l)l/2x(k+l)/2  l  r=2j  r=2T4 

j>0  j>0 


5* 


2,-i  r(k)(x-i),/Vk+,,/ 


_  r  ( k4- 1 .  1 .  Al  ( [ _ I') 

2imI  2’  2’  2  v  x7 


Ae  7  +  0  F  /k  +  1.  3. 

^rf^xk/J+1‘  ‘l2  ‘  2’  2i  xj 


=  f,(x)  -r  f:(x)  . 

fj  and  f0  are  defined  by  the  first  and  second  term  in  (23),  respectively.  It  is  easily  derived  that 


_-3  _  ~k  -f  5  Ai  U(x) 

2(x  —  1)  2x  +2x2f 


since  using  (17)  and  the  same  technique  as  that  U9ed  in  Sections  2  and  3  we  can  show  that  both  f,(x) 
and  f2( x )  satisfy  the  differential  equation  (25).  Comparing  (25)  with  (2)  we  have  m  =  2,  £j  =  1  and  i2 


2  (1  / 

=  2.  As  in  the  case  of  the  noncentral  \  (or  F),  a  recursive  formula  for  f  (x)  (i  =  2,  3,  .  .  .  )  in  terms 
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of  lower  order  derivatives  can  be  easily  obtained  from  (25).  Thus  we  need  only  the  evaluation  of  f(x) 
and  f'(x)  to  calculate  all  required  derivatives.  It  is  easily  obtained  from  (24)  that 


(V)  = 


2(7^)-^  «*>  + 


A  /2  r(k+i)  A2(k+!)  /k+3  3  a2  \ 

*  V  ’  *  *>) 


4>n?  r(j)(x-i) 7  x 


k+o  A/2r(|  +  i)(k  +  2)  /  2  \ 

2x  !<x>+  6^;  nj)  *k/!+3  1  v  +  2:  5;  2  ('-*))' 


It  follows  from  (23)  that  evaluation  of  hypergeometric  functions  of  type  jF^a;  b;  z)  is  a  main  step 
which  has  I  wen  addressed  in  Section  3. 

(2) 

Some  numerical  results  are  given  in  Table  3  to  show  the  quick  convergence  of  the  Gn  -transform. 
Like  in  the  previous  cases,  n  =  2  or  3  is  usually  sufficiently  accurate. 


5.  Conclusion 

( m ) 

In  this  pa[wr.  we  have  considered  the  applications  of  the  general  method  of  Gn  -transform  for 

tail  probabilities  to  the  three  most  commonly  used  noncentral  distributions.  The  validity  of  such 

applications  was  verified  and  accurate  numerical  results  were  given. 

Except  for  calculating  the  density  and  its  derivatives,  the  G^  ^-transform  does  not  use  any 

particular  properties  of  a  specified  distribution,  making  it  a  very  general  method.  The  method  is  also 

(rn) 

?<  ily  implemented  in  practice.  A  short  self-explanatory  FORTRAN  subroutine  for  the  Gn 
transform  is  given  in  the  appendix. 
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APPENDIX:  A  SUBROUTINE  FOR  THE  G(m> -TRANSFORM 


SUBROUTINE  GMNTRANS  (M,  N,  L,  X,  FDVEC,  GMN) 

C 

C  THIS  SUBROUTINE  CALCULATES  THE  GMN-TRANSFORM. 

C  INPUT:  M,  N,  Lf  X,  FDVEC;  OUTPUT:  GMN ( 1 )  ,  .  .  . , GMN ( N ) . 

C  M  IS  DEFINED  IN  EQUATION  (2); 

C  N  IS  THE  MAXIMUM  ORDER  OF  THE  TRANSFORM? 

C  L  -  (L(l)  ,  .  .  .  ,L(M)  )  IS  DEFINED  IN  EQUATION  (2)  (SET 

C  L(I) "I,  IF  L(I)  UNKNOWN); 

C  FDVEC  m  (FDVEC ( 1)  , FDVEC (M*N+M) )  ,  WHERE  FDVEC(I) 

C  IS  THE  (I-l)TH  DERIVATIVE  OF  THE  DENSITY  AT  X. 

C 

IMPLICIT  REAL* 8  (A-H,0-Z) 

DIMENSION  FDVEC (31)  ,  L(3)  ,  Al(30,30),  A2(30,30),  GMN(15) 
DIMENSION  WXA( 30) ,  JVEC(30,30),  PDM(30,31),  B(l) 

I JOB  *  4 
DO  100  NI=1 ,  N 
NTL  =  1+M*NI 
A1  ( 1 , 1 )  =  1 . DO 
A2 (1, 1)  =  0 .  DO 
DO  10  J  =  2,  NTL 
A1 ( 1 , J)  =  0 . DO 

10  A2 ( 1 , J)  =  FDVEC (J-l) 

DO  40  J  =  1,  NTL 
JVEC ( 1 , J )  =  1 
JVEC ( J , J )  =  1 

IF  (J.EQ. 1. OR. J.EQ. 2)  GO  TO  40 
IF  (J.EQ. 3)  GO  TO  30 
JMD2  =  J/2-1 
DO  20  X  =  1,  JMD2 

JVEC ( K+ 1 , J )  =  JVEC (K , J ) * ( J-K) /K 
20  JVEC ( J-K, J)  =  JVEC (K+ 1 , J ) 

30  JD  =  (J-l)/2 

IF  (J-1.EQ.JD*2)  JVEC ( JD+ 1 ,  J)  =  JVEC  ( JD,  J)  *  ( JD+ 1 )  / JD 
40  CONTINUE 

DO  60  K  =  0,  M-l 
DO  60  I  *  1,  NI 

POW  *  L(K+1) -1+1 

PDM (I+K*NI+1,1)  ■  X**POW 

DO  60  J  -  1,  NTL 

PDM  (I+K*NI  +  1,J+1)  -  PDM  (I  +  K*NI  +  1,J)  *(  POW-J+1 )  /X 
A1  ( I+K*NI+1 , J)  -  0 . DO 
DO  50  II  -  1,  J 

50  A1  (I+K*NI  +  1 , J)  -  A1(I+K*NI  +  1,J) + JVEC ( 1 1 , J) 

$  *PDM (I+K*N  1  +  1,11) *FDVEC ( J-1 1  +  K+l ) 

A2 ( I+K*NI+1 , J)  -  A1(I  +  K*NI+1,J) 

60  CONTINUE 

D1  -  10. DO 

CALL  LINV3F (A1 ,  B,  IJOB,  NTL,  30,  Dl,  D2 ,  WKA,  IER) 

HERE  WE  USE  IMSL  SUBROUTINE  LINV3F  TO  CALCULATE 
THE  DETERMINANT  OF  A  MATRIX. 

OTHER  EFFICIENT  SUBROUTINES  CAN  ALSO  BE  USED  HERE. 

D3  -  10. DO 

CALL  LINV3F ( A2 ,  B,  IJOB,  NTL,  30,  D3 ,  D4 ,  WKA,  IER) 

100  GMN (NI)  -  D3/D1*2 .DO** (D4-D2 ) 

RETURN 

END 
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Table  1.  Relative  error  =  R[approx.]  =  {  |  approx. —true  l/true} 


of  approximations  to  the  upper  tail  of 


k 

A 

X 

true  value 

R[gSj)] 

R[G(22)] 

R[g(33)] 

5 

1 

8 

.2466 

1.1(2)* 

1.5(4) 

1.1(4) 

13 

.07332 

2.8(3) 

1.5(5) 

4.5(7) 

21 

.003156 

4.6(4) 

5.1(7) 

3.8(9) 

5 

10 

20 

.2189 

6.2(2) 

7.9(4) 

8.3(6) 

26 

.07497 

1.9(2) 

1.5(4) 

8.4(7) 

43 

.001521 

2.3(3) 

5.4(6) 

5.5(9) 

25 

10 

44 

.1674 

8.1(2) 

8.0(3) 

6.6(4) 

50 

.06824 

3.2(2) 

1.8(3) 

9.9(5) 

70 

.001337 

3.7(3) 

4.8(5) 

7.8(7) 

25 

25 

62 

.1616 

9.5(2) 

9.4(3) 

7-4(4) 

70 

.06159 

3.7(2) 

2.0(3) 

1.0(4) 

94 

.001295 

4.7(3) 

6.4(5) 

1.1(6) 

*  1.1(2)  =  1.1 

x  10~2 
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Table  2.  Relative  error  =  R[approx.]  =  {  |  approx.  —  true  |  /true} 
of  approximations  to  the  upper  tail  of  Ff  (A). 


k 

P 

A 

y 

true  value 

R[G(2j)] 

R[g(32)] 

Edgeworth 

2 

4 

2 

5 

.2003 

2.6(4)* 

1.4(7) 

1.7(5) 

10 

.07906 

2.1(5) 

4.2(9) 

2.9(3) 

33 

.01065 

2.1(7) 

5.1(12) 

6.1(3) 

2 

4 

15 

20 

.2135 

9.1(4) 

4.9(7) 

4.2(4) 

42 

.06946 

3.3(5) 

5.5(9) 

8.1(3) 

106 

.01353 

6.7(7) 

2.3(11) 

7.6(3) 

2 

20 

10 

10 

.1817 

1.5(2) 

1.3(5) 

1.5(3) 

15 

.05327 

1.7(3) 

6.4(7) 

3.0(3) 

25 

.005622 

1-2(4) 

5.3(8) 

5.9(3) 

20 

20 

20 

3 

.1775 

1-7(2) 

8-1(4) 

4.5(4) 

4 

.05722 

2.2(3) 

5.3(5) 

6.6(5) 

10 

.002323 

1.0(5) 

3.0(8) 

12(2) 

*  2.6(4)  =  2.6  x  10-4 
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Table  3.  Relative  error  =  R[approx.]  =  {  |  approx.  — true  |/true} 
of  approximations  to  the  upper  tail  of  tj^A). 


It 

A 

y 

true  value 

r[gS2)] 

r[g(2j)] 

R[g(3j)] 

3 

1 

2 

.2564 

2.5(2)* 

1.1(3) 

1.2(4) 

4 

.06390 

1.3(3) 

3.8(6) 

5.7(8) 

8 

.01012 

7.3(5) 

1.0(8) 

1.9(11) 

3 

7 

12 

.2087 

2.9(2) 

5.4(4) 

7.1(4) 

19 

.06397 

3.4(3) 

3.6(5) 

1.9(5) 

37 

.009575 

2-0(4) 

1.4(8) 

1.9(9) 

10 

7 

10 

.1263 

5.2(2) 

1-7(3) 

3.5(5) 

12 

.04187 

1.6(2) 

2.2(4) 

2.1(6) 

17 

.003113 

2.5(3) 

6.9(6) 

1-3(8) 

*  2.5(2)  =  2.5  x  10“3 
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